Non-lattice simulation for supersymmetric gauge theories in one dimension 
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Lattice simulation of supersymmetric gauge theories is not straightforward. In some cases the 
lack of manifest supersymmetry just necessitates cumbersome fine-tuning, but in the worse cases the 
chiral and/or Majorana nature of fermions makes it difficult to even formulate an appropriate lattice 
theory. We propose to circumvent all these problems inherent in the lattice approach by adopting 
a non-lattice approach in the case of one-dimensional supersymmetric gauge theories, which are 
important in the string/M theory context. 



PACS numbers: ll.lO.Kk; 11.15.Ha; ll.30.Pb 

Introduction. — Lattice gauge theory, together with 
the developments of various simulation techniques, has 
provided us with a powerful non-perturbative method to 
study gauge theories such as QCD. However, when one 
tries to apply the method to supersymmetric gauge the- 
ories, which are interesting for many reasons, one has to 
face with some practical and theoretical obstacles. 

First of all, since the algebra of supersymmetry con- 
tains continuous translations, which are broken to dis- 
crete ones, it seems unavoidable to break it on the lattice. 
Then, one has to include all the relevant terms allowed 
by symmetries preserved on the lattice, and fine-tune the 
coupling constants to arrive at the desired supersymmet- 
ric fixed point in the continuum limit. Recent progress 
(See ref. [l|] and references therein.) is that one can re- 
duce the number of parameters to be fine-tuned (even to 
zero in some cases) by preserving some part of supersym- 
metry. In lower dimensions, one can alternatively take 
the advantage of super-renormalizability, and determine 
the appropriate counter-terms by perturbative calcula- 
tions to avoid fine-tuning. These two approaches can be 
nicely illustrated in one dimension by the example of a 
supersymmetric anharmonic oscillator [sl . 

In the string/M theory context, supersymmetric gauge 
theories appear in various ways. In particular, the low- 
energy effective theory of a stack of N p-branes are given 
by {p + l)-dimensional U(iV) super Yang-Mills theory. 
This led to various interesting conjectures. For instance, 
the gauge/gravity duality states that the strong coupling 
limit of large- TV gauge theories has a dual description in 
terms of classical supergravity. A different but closely 
related conjecture asserts that non-perturbative formu- 
lations of superstring/M theory can be given by matrix 
models, which can be obtained by dimensionally reducing 
lOd A/" = 1 U(iV) super Yang-Mills theory to D = 0, 1, 2 
dimensions. In particular, the D = 1 model [3| corre- 
sponds to the M Theory. Its bosonic version has been 
studied by Monte Carlo simulation using the lattice for- 



mulation [5| and the continuum quenched Eguchi-Kawai 
model 0. However, the next step of including fermionic 
matrices would be difficult even theoretically because of 
their Majorana- Weyl nature inherited from lOd. See ref. 
\1] for a proposal using lattice, which preserves half of 
SUSY at the expense of breaking the SO (9) symmetry. 

The aim of this letter is to point out that there ex- 
ists an extremely simple and elegant method to simulate 
supersymmetric gauge theories in one dimension. (We 
discuss the case of U(A^) gauge group, but extension to 
more general group must be possible.) In particular, it 
will enable us to study the gauge theory side of the afore- 
mentioned duality and also to put M theory on computer. 

Let us first recall that the importance of the lattice 
formulation lies in its manifest gauge invariance. In the 
present Id case, however, the gauge dynamics is almost 
trivial. (We assume that the Id direction is compact. 
The non-compact case would be easier since the gauge 
dynamics is completely trivial.) This gives us an oppor- 
tunity to use a non-lattice formulation. More specifically, 
we first take the static diagonal gauge. Using the residual 
large gauge transformation, we can choose a gauge slice 
such that the diagonal elements of the constant gauge 
field lie within a minimum interval. Finally we expand 
the fields into Fourier modes, and keep only the modes 
up to some cutoff. The crucial point of our method is 
that the gauge symmetry is completely fixed (up to the 
global permutation group, which is kept intact) before 
introducing the cutoff. This is specific to one dimen- 
sion, and the momentum cutoff regularization in higher 
dimensions generally breaks gauge invariance. 

Supersymmetric anharmonic oscillator. — To gain 
some insight into our new approach, we first apply it to a 
non-gauge supersymmetric theory, which is well studied 
by the lattice formulation. In particular, supersymmetry, 
which is broken by the cutoff in our formalism, is restored 
much faster than the continuum limit is achieved. 

While the manuscript was being prepared, we received 
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a preprint [g], in which the same model is studied on 
the lattice using various methods. As far as non-gauge 
theories are concerned, our approach is almost equivalent 
to the method with the non-local SLAC derivative 
The only difference is the identification of the modes at 
the boundary of the Brillouin zone in the lattice case. As 
a consequence, our results shown in fig. [1] agree with the 
corresponding results in ref . Q . 
The model is defined by the action 
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dt 



{{dtcl>f + h'{cbf}+i^{dt + h"{cl,))ij 



(1) 

where (/> is a real scalar field, and ip is a one-component 
Dirac field, both in Id, obeying periodic boundary condi- 
tions. This model is supersymmetric for arbitrary func- 
tion h{(j)), but here we take = ^m(jp + \g<l>^- In our 
approach we make a Fourier expansion 
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and similarly for the fermionic fields, where n takes inte- 
ger values, and A is the UV cutoff. In terms of the Fourier 
modes, the action can be written as S* = S'b + S-p, where 
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We have introduced a shorthand notation 
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Integrating out the fermions first, we obtain the effective 
action for the bosons as 



S^ff = 5b - IndetX 



(6) 



where the fermion determinant dct M. is real positive for 
positive m and g. 

As an efficient algorithm to simulate the model (O, 
we use the idea of the hybrid Monte Carlo algorithm 
[lo| . Let us introduce the auxiliary real field n(<), whose 
Fourier components are denoted as n„ 
action 



and consider the 



5, 
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Integrating out the auxiliary field first, we retrieve 
In order to update the fields, we solve the equations 



dT 
dT 



Q!„— = anli-n 



(8) 
(9) 



along the fictitious time r for a fixed interval Tf . The real 
coefficients q;„ should be optimized based on the idea of 
the Fourier acceleration This evolution, if treated 

exactly, conserves the action S'hmc- In practice, we dis- 
cretize the r-evolution in such a way (the leap-frog dis- 
cretization) that the reversibility is maintained. Due to 
the discretization, the action 5hmc changes by a small 
amount, say AS'hmc- In order to satisfy the detailed 
balance, we accept the new configuration with the prob- 
ability min(l, e"'^'^'^"'^), which is the usual Metropolis 
procedure. Before we start a new r-evolution, we refresh 
the n„ variables by drawing Gaussian random numbers 
which follow from the action ([7]). This procedure is nec- 
essary for avoiding the ergodicity problem. (The step size 
Ar should be optimized for fixed rf by maximizing the 
acceptance rate times At. Then n should be optimized 
by minimizing the autocorrelation time in units of step 
in the r-evolution.) 

The main part of the computation is the evaluation of 
the term in eq. ([9]) given by 
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The convolution requires O(A^) calculations, while the 
inverse M^^ requires ©(A"^) calculations. 

As usual, we extract masses from the exponential decay 
of the two-point functions 

A 

GB{t) = ((^(0)0(0) = 6o + 2 ^ 6„ cos(wni) , (12) 

n=l 

A 
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where we have defined 6„ = {\4>n\^) and c„ = {{M~^)nn) ■ 
For the fermion, it proved convenient to consider , instead 
of (fT3|) . a symmetrized correlator 

4^^"^)(t) ^ }-{Gp{t) + GA-t)} 



Co 



2 ^ Re(c„) cos(wni) , (14) 



where we have used the fact {M.nk)* = M.-n.-k- In fact 
the functions and with respect to t oscillate with 



3 



the frequency of the order of cutoff. This is nothing but 
the Gibbs phenomenon due to the sharp cutoff in the sum 
over Fourier modes. To overcome this problem, we note 
that the coefficients &„ behave as 5„ ~ ^ + at large n 
as can be shown from general arguments. We obtain the 
coefficients di and d2 from the results at n = A — f, A, 
and extend the sum in (|f 2p over n up to f 000 assuming 
the above asymptotic form. We make an analogous anal- 
ysis for Re(cn) in (fT4|) . In this way we are able to see 
clear exponential behaviors, and extract the correspond- 
ing masses. The results for /3 = 1, m = 10, g = 100 
are plotted against 1/A in fig.[TJ (Note that the effective 
coupling constant is g/m^ = 1-) We find that the finite 
A effects are 0(1/A), and that the data points for the 
boson and the fermion lie on top of each other. Thus 
in our formalism, the effect of supersymmetry breaking 
by the cutoff disappears much faster than 1/A. In the 
same figure we also plot the results obtained from lattice 
formulations for comparison. (Matching the number of 
degrees of freedom, we make an identification A = ^, 
where L is the number of sites.) 




FIG. 1: The circles (diamonds) are the mass for the boson 
(fermion) obtained by our method for A = 8, 10, 12, 14, 16. 
The triangles (inverted triangles) are the mass for the boson 
(fermion) obtained by Giedt et al. [Jl with the 0(a)-improved 
lattice action, and the squares are the results obtained by Cat- 
terall and Gregory with the lattice action preserving half 
of SUSY, hence degenerate. The horizontal line represents 
the exact result, and the dotted lines represent the expected 
behaviors at large A. 

Supersymmetric matrix quantum mechanics. — Here 
we consider a model with four supercharges defined by 



S 



1 



[Xi, Xj 



(15) 



where Dt — dt ~ i[A{t), ■ ] represents the covariant 
derivative with the gauge field A{t) being a,n N x N Her- 
mitian matrix. The bosonic matrices Xi{t) (i = 1,2,3) 



are N x N Hermitian, and the fermionic matrices ipa{t) 
and -ipait) (a = 1,2) are N x N matrices with complex 
Grassmann entries. The 2x2 matrices ct^ are the Pauli 
matrices. The model can be obtained formally by di- 
mensionally reducing 4d A/" = 1 U(iV) super Yang-Mills 
theory to Id, and it can be viewed as a one-dimensional 
U(A^) gauge theory. (The totally reduced model has been 
studied by Monte Carlo simulation in refs. [l2|-) Let us 
assume the boundary conditions to be periodic for bosons 
and anti-periodic for fermions. The extent (3 in the Eu- 
clidean time direction then corresponds to the inverse 
temperature /? = 1/T. The parameter g in (jl5|) can al- 
ways be absorbed by an appropriate rescaling of the ma- 
trices and the time coordinate t. 
without loss of generality. 

Let us take the static diagonal gauge A{t) = 
■|diag(ai, • • ■ , ctAr), where aa {a — 1, ■ ■ ■ , N) can be cho- 
sen to lie within the interval (— 7r,7r] by making a g auge 
transformation with a non-zero winding number [l3| . We 
have to add to the action a term 



Hence we set g = 
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which appears from the Faddeev-Popov procedure, and 
the integration measure for aa is taken to be uniform. 
We make a Fourier expansion 
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and similarly for tp, where r takes half-integer values, 
due to the anti-periodic boundary conditions, and A = 
A — 1/2. Eq. (fTS]) can then be written as 
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The algorithm for simulating ([T8|) is analogous to the 
previous model. Here we introduce the auxiliary vari- 
ables Ili{t) and Pa, which are N x N Hermitian matri- 
ces conjugate to Xi{t) and N real variables conjugate to 
aa, respectively. The fermion determinant is real pos- 
itive, and the computational effort for one step in the 
r-evolution is proportional to A'^A^^. Figures [2] and [3] 
show the results for the energy and the Polyakov line, 
respectively, for the bosonic and supersymmetric cases. 
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FIG. 2: The energy (normalized by A^^) is plotted against 
temperature for the matrix quantum mechanics with A'^ = 4. 



be quite useful. In particular, the Fourier acceleration 
requires no extra cost, since we deal with the Fourier 
modes directly. The continuum limit is achieved much 
faster than one would expect naively from the number 
of degrees of freedom. This is understandable since the 
Fourier modes omitted by the cutoff scheme are naturally 
suppressed by the kinetic term. 

It is straightforward to apply our method to the most 
interesting case with sixteen supercharges, which is cur- 
rently under investigation [ist . 
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FIG. 3: The same as fig. [2] but for the absolute value of the 
Polyakov line. 

In the bosonic case we also plot the results from lat- 
tice simulation with the lattice spacing a = 0.02. (The 
number of lattice sites is given by L = 1/ (Ta) , which is 
50 for T = 1.) The results obtained by our new method 
approach the lattice result as A is increased. 

In the supersymmetric case, our preliminary results 
with A 8 reproduce the asymptotic behavior at large T 
obtained by the high temperature expansion (HTE) 
up to the next-leading order. (The solid lines represent 
the results at the leading order of HTE, which are the 
same for the bosonic and SUSY cases.) Note that our 
method is applicable also at low temperature, where the 
HTE is no more valid. 

Summary and concluding remarks. — In this letter we 
have proposed a new simulation method, which enables 
non-perturbative studies of supersymmetric gauge the- 
ories in one dimension. For practical implementation, 
the idea of the hybrid Monte Carlo algorithm seems to 
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